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Abstract. We perform a detailed quantitative analysis of the recent AMS-02 electron and 
positron data. We investigate the interplay between the emission from primary astrophysical 
sources, namely Supernova Remnants and Pulsar Wind Nebulae, and the contribution from 
a dark matter annihilation or decay signal. Our aim is to assess the information that can 
be derived on dark matter properties when both dark matter and primary astrophysical 
sources are assumed to jointly contribute to the leptonic observables measured by the AMS- 
02 experiment. We investigate both the possibility to set robust constraints on the dark 
matter annihilation/decay rate and the possibility to look for dark matter signals within 
realistic models that take into account the full complexity of the astrophysical background. 
Our results show that AMS-02 data enable to probe efficiently vast regions of the dark matter 
parameter space and, in some cases, to set constraints on the dark matter annihilation/decay 
rate that are comparable or even stronger than the ones derived from other indirect detection 
channels. 
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1 Introduction 

Dark matter (DM) stands out as one of the most intriguing challenges in physics. While 
gravitational effects of this elusive substance are manifest in a variety of observables, ranging 
from the astrophysical to the cosmological scale, an indisputable non-gravitational evidence 
of its existence has yet to be unveiled. Among the techniques developed to pursue the search 
of DM non-gravitational imprints, there is the so-called indirect detection, which consists in 
the attempt to reveal a signal that can be the outcome of a DM pair annihilation or decay 
event and that can possibly manifest its presence as an exotic component in cosmic rays 
(CRs). 

In the framework of this searching strategy, a preeminent role is played by the analysis 
of the fraction of antimatter in the flux of charged CRs. In particular, in recent years, the 
observations performed by the PAMELA satellite of a steep rise in the energy spectrum of 
the positron fraction [1] has been interpreted by many authors as a possible evidence for 
the existence of a leptophilic TeV-scale WIMP (see, for example, Refs. [2-9]). However, as 
stressed in several works, plausible astrophysical mechanisms able to generate such a rise in 
the positron fraction can be conceived: while the most widely known example invokes Pulsar 
Wind Nebulae (PWNe) or Supernova Remnants (SNRs) or a combination of both as sources 
of primary [10-14], it has also been shown that the secondary production of within the 
shock region of SNRs provides a perfectly viable mechanism to explain observations [15, 16]. 
Several features that characterize these interpretations of the PAMELA excess have been 
investigated in a series of recent works (see, for example. Refs. [17-27]) in connection to the 
data delivered by the AMS-02 experiment. 
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In a previous work [28] we have carried out a detailed and quantitative analysis of the 
four observables measured by the AMS-02 experiment (electron flux, positron flux, total lep- 
tonic flux and positron fraction) within a purely astrophysical model for generation. In 
this paper our purpose is to extend this quantitative study to include also a DM contribu¬ 
tion: the aim is to assess the constraints on DM obtainable within a model in which the 
astrophysical background is included and modeled in a realistic way. In particular, we stress 
that, differently from Ref. [28] where we used an 8% error on the AMS-02 measured fluxes, 
we are now making full profit of the recently published AMS-02 Collaboration data, which 
report smaller uncertainties on their whole set of observables. 

The paper is organised as follows: in Section 2 we briefly illustrate the various contri¬ 
butions to the electron and positron fluxes, while in Section 3 we describe the framework 
that we use to model the transport of positrons and electrons across the Galaxy and the 
heliosphere. In Section 4, we present the methods that we use for the different parts of our 
analysis and the results that we obtain. Finally, in Section 5 we derive our conclusions. 


2 Contributions to the e'''e flux 


In this Section, we briefly review the role of the main sources that contribute to the four 
observables measured by the AMS-02 experiment. For a more detailed description of all the 
astrophysical contributions (primary and secondary) we address the reader to Ref. [28]. 


2.1 Secondaries 

Electrons and positrons can be generated in spallation reactions that involve primary CRs, in 
particular proton and Helium nuclei, impinging on the nuclei that populate the interstellar 
medium (ISM), the dominant component of which being represented by Hydrogen. The 
source term related to this contribution is: 


q^±{x,Ee) = dvr nisM(^) 


dEcR^CR {x, Ecr) -^{Ecr, Ee), 


( 2 . 1 ) 


where the quantity nisM represents the density of target nuclei in the ISM, <I>cr is the flux 
of the primary CR species and the term da/dEf, stands for the differential inclusive cross 
section for the electron/positron production in the spallation reaction under study. 

By following the same approach outlined in Ref. [28], we set the primary CR fluxes by 
fitting the recent AMS-02 measurements of the proton [29] and Helium [30] fluxes with the 
following interstellar spectra: 


‘I’cr(-R) = 


for R < Rbreak, 

for R > Rbreak, 


( 2 . 2 ) 


where R = pc/eZ and (3 denote, respectively, the rigidity and the velocity of the particle 
under consideration. We allow for the presence of a spectral break at a rigidity Rbreak 
to be determined through the fitting procedure. For the best-fit configuration the values 
of the parameters are: A = 26700 ± 500 m“^s“^sr“^(GeV/n)“r^ p = 7.2 ± 0.4, Pi = 
2.877±0.004, P 2 = 2.748±0.013 and Pbreak = 220±22 GV for the proton and A = 4110±80 
m-2s-rsr-r(GeV/n)-i, P = 3.5 ± 0.7, Pi = 2.793 ± 0.004, P 2 = 2.689 ± 0.013 and Pbreak = 
187 ± 18 GV for the Helium, the Fisk solar modulation potential being set to 700 MV. The 
best fit to the combined proton and Helium data has a reduced chi-square (x^/d.o.f.) equal 
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Fit to AMS-02 p and He data 



Figure 1. Fit to the proton and helium fluxes measured by AMS-02: the black points denote 
experimental data, while the blue and red solid lines represent our best fits. 


to 0.43.The proton and Helium fluxes for this best-fit configuration are shown, together with 
the AMS-02 measurements in Fig. 1. 

Concerning the inclusive production cross sections da/dEf., for the case of p-p collisions 
we use the parameterization described in Ref. [31], while to model processes involving Helium 
nuclei, either as the incoming particle or the target, we resort to the empirical prescriptions 
illustrated in Ref. [32]. 


2.2 Supernova Remnants 

SNRs are usually thought to play the dominant role in accelerating charged cosmic rays 
within the Galaxy. This acceleration is realised through the propagation of non-relativistic 
shock waves that are produced by the star explosion. The typical SNR injected spectra have 
the shape of a power-law with a cut-off at large energies: 

/ E \ "'^SNR= / E \ 

Q { E ) = Qo,snrs ( I \ ~^) ’ 

being Qo.snrs the normalization of the spectrum, 7 snrs the spectral index and Ec the cut-off 
energy. As it has been extensively discussed in [28], the parameters (5o,snrs and 7snrs can 
be inferred from the study of the radio emission measured in the region of the sky that hosts 
the SNR, while E^ is expected to lie at the TeV range [33-37]. In the analysis that will 
be described in the following Sections of the paper, we assume Ec = 2 TeV. Let us remark 
that, however, as far as Ec is beyond the maximal energy that is measured by AMS-02, the 
exact value of this parameter does not affect our results. Finally, every time that, in our 
investigations, we will need to use SNRs parameters we will resort to the Green catalogue 
[38], which is the most complete catalogue of Galactic SNRs. 
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2.3 Pulsar Wind Nebulae 


Also pulsars are expected to generate a flux of electrons and positrons through a mechanism 
known as spin-down emission (for a description of which we address the reader to Refs. [39- 
47]). In a few simple words, this mechanism is a consequence of electrons and positrons being 
torn away from the surface of the neutron star by the strong electric field genarated by the 
pulsar spinning. These charged particles gather in a sort of wind that surrounds the pulsar 
and then are released in the ISM at the disruption of this nebula. Because of this injection 
mechanism, which is fast and followed by a weak residual energy emission, PWNe can be 
considered as burst-like sources of e^. 

The spectrum of the injected by a PWN in the ISM has the same expression as the 
one in Eq. 2.3 associated to SNRs. As outlined in [28], the normalization of this spectrum is 
related to the efficiency r/pwNe with which the PWN can convert its spin down energy into 
the production of pairs: 



dEEQ{E) 


?/PWNe kPO) 


(2.4) 


where the quantity Wq represents the total spin-down energy which, in terms of the present 
age of the pulsar E and the typical pulsar decay time tq can be expressed as: 


IPo 


tqE 



(2.5) 


The most complete list of PWNe is represented by the ATNF catalogue [48]. As it will 
be widely discussed in the following, we will use it as a reference for all the PWN parameters. 


2.4 Dark Matter 

Positrons and electrons can also be the result of the pair annihilation or decay of DM particles. 
The source terms associated to these contributions are: 


Qann(^, E) = 


Qdec(^) E^ — 



( 2 . 6 ) 


where x denotes Galactic position, e being a factor that takes the value 1/2 or 1/4 for, re¬ 
spectively, a self-conjugate or non self-conjugate DM particle, while / denotes the Standard 
Model particles that can be produced in the annihilation or decay process and the functions 
dN^^/clE represent the energy spectrum generated in a single annihilation or decay pro¬ 
cess. The galactic DM halo, filled with particles with mass rriDM-, follows a spatial density 
p{x). We perform a model independent analysis which consists in assuming that the DM 
annihilation/decay occurs in a single channel. In particular, we will focus our attention on 
the five channels e+e", r+r“, bb, W~^W~. We model the energy spectra dN^^/dE 

from Ref. [49]: we remind that these spectra have been computed by taking into account 
electroweak correction which, as stressed in [50] can play a non-negligible role in shaping the 
emission when the DM mass is above the electroweak scale. 
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3 Transport in the galaxy and in the Heliosphere 

After being injected by their source, electrons and positrons propagate across the ISM where 
spatial diffusion, convection, reacceleration and energy losses can shape their spectrum. To 
model this transport and take into account these effects, we have to solve the transport 
equation: 

dtM-V- {k{E)VM] + I = Q{E, X, t) , (3.1) 

where N denotes the positron(electron) number density per unit energy while Q{E, x, t) 
represents the source terms that were discussed in the previous Sections. The diffusion 
coefficient K{E), which we assume to be space independent, is defined in the following way: 

K{E) = /3Ko{R/l GYY- (3.2) 

The term ^ is given by the sum of the synchrotron energy loss term (which depends on the 
value of the Galactic magnetic field, for which, compatibly with Ref. [51], we have assumed an 
intensity of 2 /rG (3 fiG) for the regular (random) component) and the inverse Compton (IC) 
loss term (for which we have included a full relativistic treatment [52]). Additional energy 
loss processes, namely bremsstrahlung, ionization and Coulombian interactions on the ISM 
can be safely neglected. 

In solving Eq. 3.1, we adopt the two-zone diffusion model, which has been extensively 
described in literature (see for example Ref [53]). In this paper we adopt the set of propagation 
parameters defined in the Med propagation model of Ref. [54]. For any further detail we 
address the reader to Ref. [52]. 

Lastly, for the solar modulation of the fluxes, we have used the usual force field approx¬ 
imation treated in Refs. [55, 56]. In this frameworks, the solar modulation depends on the 
value of the parameter (f>, the so-called Fisk potential. As it will be discussed in the following, 
we will always determine (p directly from data. 


4 Analysis and results 


The AMS-02 Collaboration has recently published the results of their 30 months measure¬ 
ments of the full set of observables that are related to the physics of cosmic electrons and 
positrons, i.e., the e~ and e+ fluxes [57], the total e~ +e'^ flux [58] and the positron fraction 
e"*"/(e“ -l-e''') [59]. Our main purpose in this Section is to test various possible interpretations 
of these experimental data. 

The analyses are based on a global fit of the four aforementioned observables within 
several theoretical frameworks (described below), each one being characterized by a certain 
set of parameters {6i,... 6^}, whose total number N depend on the specific model under 
scrutiny. The iV-dimensional parameters space of each model is sampled by means of a 
Markov Chain Monte Carlo (MCMC) scan performed with the COSMOMC package [60]. For 
every parameter considered in our analyses, a flat prior is assumed. In order to determine 
the configuration of parameters providing the best agreement with AMS-02 measurements 
we assign to each N-tuple {9i,... 9 m} a chi-square estimator defined as: 


X 


= — 21nL = 


EE 


ifi i9i, ■ ■ ■ 9n) - dlY^ 


■J\2 


(4.1) 



where the indices j and i run over the four AMS-02 observables and the different energy 
bins, respectively. The quantities fl{0i,.. .On) represent the theoretical predictions while 
dl stands for the observed experimental data values, which are affected by the experimental 
uncertainties 

The experimental uncertainties al are defined as the sum in quadrature of the statistical 
and systematic uncertainties. The latter, reported in [57-59], are associated to event selection, 
determination of the acceptance, bin-to-bin migration and charge confusion (that, obviously, 
affects all the observables except the total e"*“ + e~ flux) and in the interval of energies that 
we consider in our analysis have an impact that ranges from a few % (for energies around 10 
GeV) to 10-15 % (for the highest energies). 

Another point worth to be mentioned here is that in this analysis we choose to consider 
only data above 10 GeV. In fact, the sector relevant for the investigation of primary astro- 
physical sources and of the corresponding DM contribution is the high energy part of 
the spectra. Gonsidering that time-dependent solar modulation effects become increasingly 
important at lower energies, we decide to concentrate our analysis only on the high energy 
part, in order to avoid that the numerous AMS-02 energy bins at very low energy (which 
require detailed, although uncertain, solar modulation modeling) would drive the statistical 
significance of the fit to the part of the spectrum which is of least importance to the present 
analysis. Since the effect of transport in the heliosphere is still present up to few tens of GeV, 
we nevertheless consider solar modulation by applying the simple Fisk method (as discussed 
above): however, by concentrating on energies above 10 GeV we minimize the bias of the 
low-energy solar-modulation-dominated data bins. The choice of 10 GeV as the discriminat¬ 
ing energy has been adopted as a trade off between the requirement of reducing the impact 
of solar modulation modeling and the interest to investigate the energy region where DM 
and primary sources start to emerge. 

4.1 The purely astrophysical interpretation 

In this paragraph, we follow the prescriptions presented in our previous work [28] and discuss 
our new analysis for a purely astrophysical interpretation of AMS-02 measurements. In brief, 
the model, which we label astro model, is composed by: 

• Secondary e~ and e"*": their flux is computed by using the source term introduced in 
Section 2.1. We recall that this component is fixed. In fact, its contribution depends 
on the primary GR fluxes and on the parameterization of the spallation cross sections, 
which we do not vary in our analysis. The secondary e"*" components is the dominant 
one at low energies, where it has been shown in Ref. [28] to perfectly reproduce the 
measured AMS-02 flux, within the chosen propagation model. 

• Primary e~ accelerated by SNRs: for this contribution, which has been discussed 
in Section 2.2, we follow the approach outlined in Ref. [28], that consists in separating 
the total population of SNRs into two categories: 

— Far SNRs {d > 3 kpc): they are modeled as an average component with a spatial 
distribution taken from Ref. [61]. The sources that fall in this class share common 
values for the normalization (5 o,snrs and the spectral index 7 snrs of the injected 
electron flux. These parameters are left free to vary in our fitting procedure. 

— Near SNRs (d < 3 kpc): the parameters that characterize these sources, i.e. their 
distances, ages, spectral indices and the normalization of their fluxes, are taken 
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directly from the Green catalogue [38], as anticipated in Section 2.2 and extensively 
discussed in Ref. [28]. Notice that, in order to allow for a better agreement with 
the electron flux at high energies, we introduce a free normalization for the Vela 
SNR, Vveiaj which therefore represents an additional parameter, not included in 
our previous analysis of Ref. [28]^. Freedom in the variation of this parameter 
can be ascribed to variations in the value of the magnetic field of the source: 
Vveia oc R“(t'snrs+i)/ 2^ Vyeia = 1 Corresponding to R = 30 /iG. 

• Primary e~ and e"*" accelerated by PWNe: by using the same approach of Ref. [28], 
we consider the PWNe of the ATNF calatog [48] as primary sources. The pulsars spin- 
down energies are derived from the catalog itself. The efficiency ?/pwNe of conversion of 
this energy into electrons and positrons, and the spectral index ypwNe of the emission 
are assumed to be common to all the pulsars in the catalog and are free to vary. 

The contributions to the fluxes coming from all these sources are propagated in the 
ISM and modulated in the heliosphere according to the prescriptions discussed in Section 
3. While the reference model for the Galactic propagation is described by the Med set of 
parameters, the Fisk potential cj) of the solar modulation is a free parameter. To summarize, 
the astro model is characterized by six free parameters: Qo.SNRs) Tsnrsj -Vveiaj f/PWNe, 7PWNe 
and 4>. 

We start our analysis of the astro model by fitting the whole amount of data released by 
AMS, which consists of the four datasets corresponding to the positron flux, the electron flux, 
the total e"*" + e~ flux and the positron fraction. Fig. 2 illustrates the posterior probability 
distributions of the parameters of the model, as they are sampled by the MGMG scan. Their 
values for the best-fit configuration are reported in the first column of Table 1 , together with 
the global chi-square and the chi-square associated to each dataset. While the overall fit is 
good, the contribution to the total chi-square that comes from the positron fraction is large 
(the per data point is close to 2). This points toward a hint that the astro model might 
have some issues in reproducing the behaviour of this observable. 

To have a clearer view about this fact, we separate the four datasets in two independent 
groups, and fit them separately: in one case we use the single electron and positron fluxes, 
in the other case the positron fraction and the total electron-l-positron flux. The second 
column of Table 1 reports the best-fit parameters obtained by fitting only the positron frac¬ 
tion and the total flux, while the third column illustrates the outcome of a fit performed 
only on the positron and electron fluxes. The best-fit configurations for the former case are 
shown, together with AMS-02 data, in Fig. 3. The results confirm that the astro model 
has no problem in predicting single and total fluxes that are in a agreement with the ones 
measured by AMS-02, but it has some tension with the measured positron fraction^. Notice 
that the AMS-02 data are internally perfectly consistent, as can be seen in Fig. 4 where the 
positron fraction and total flux derived from the single fluxes are plotted together with the 
corresponding measured quantities. We have explicitly checked that the deviation between 
the reconstructed and measured quantities is always well within the uncertainty of the recon¬ 
structed values, which is determined by adopting the usual error propagation prescriptions 
(the only exception being the total flux at 262 GeV which however does not represent an 
issue at all in our analysis). The main difference arises from the much smaller experimental 

^We keep the distance and age of Vela fixed to, respectively, dveia = 0.294 kpc and Tveia = 11-3 kyr 
^We have checked that analogous results are obtained if different combinations of datasets are considered. 
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Figure 2. Triangular plot for the fit to the four AMS-02 datasets of the astro model parameters, as 
reported in the first column of Table 1. The contours refer to Icr, 2cr and Scr C.L. allowed regions. 
The plots along the diagonal show the posterior distribution for each parameter. 


uncertainty on the directly-measured position fraction (especially for energies below 30 GeV), 
as a consequence of systematic effects that cancel out in the determination of ratios. 

This fact is further investigated in Fig. 5: in the top left panel we plot the relative 
uncertainty associated with the different datasets as a function of the energy. The curve for 
the positron fraction is steeper than the other curves and, for energies below 30 GeV it is up 
to a factor of 3 lower than the uncertainties characterizing the other datasets. Notice also 
that the relative uncertainty for the positron fraction is always smaller than the one for the 
positron flux, in the whole energy range. Since positrons are the more relevant observables 
to study dark matter, and since the positron fraction is experimentally better determined, 
we decide to adopt this as the key observables for our analysis. And since this observables 
appears to have some tension with our pure astro model, we investigate how the agreement 
can be improved. 

In this regard, it is useful to recall that, within the assumptions that characterize the 
astro model, the behaviour of the positron fraction is determined by the two parameters 
related to the emission from PWNe. While the effect of ??pwNe is basically a normalization, 
a change in the ypwNe parameter reflects into a change in the slope, as shown in the top 
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0.76 

xlf (43 data pts) 

81.7 

80.4 

- 

xlum (50 data pts) 

36.0 

37.2 

- 

Xg+ (49 data pts) 

39.7 

- 

40.4 

Xg_ (49 data pts) 

33.6 

- 

28.9 


Table 1. Best-fit parameters for the astro models together with the associated to the different 
datasets, for different combinations of the datasets. 



E [GeV] 


e*/(e*+e) 



E [GeV] 


Figure 3. Astrophysical fit {astro model) to the AMS-02 positron fraction and total flux datasets 
[57, 58, 62]. The fit refers to the second column of Table 1. We display also Fermi-LAT [63, 64], 
PAMELA [1, 65, 66], and HESS data [33, 67]. The styles and colors used to represent the various 
contributions are described in the insets. 


right panel of Fig. 5. It is quite evident that deviations of qpwNe from the best-fit value 
1.95 are strongly constrained by the large number of data points between 20 and 100 GeV. 
And obviously, once that ypwNe is set, fluctuations in the normalization r/pwNe away from 
the best fit value 0.037 are not possible. 

To illustrate how small variations in these two parameters affect the quality of the fit, we 
plot in the bottom panel of Fig. 5, for each data point of the positron fraction, the chi-square 
that is associated with the best-fit configuration of the astro model (for the combined fit of 
the positron fraction and of the total flux) together with the chi-square for slightly different 
values of the parameters of the PWNe. As it can be seen, the data points giving the maximal 
contribution to the total chi-square are the two at the lowest energies and the two at 77.1 
GeV and 123.31 GeV (this arises from the wiggly features that the positron fraction possesses 
at these energies). It is evident that, in the astro model, every change of parameters that 
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Figure 4. Comparison between the directly measured total flux and positron fraction (black points) 
and the corresponding observables reconstructed from the measured single positron and electron fluxes 
(red points). 


goes in the direction of lowering the chi-square at low energies has the effect of augmenting 
the chi-square at high energies. This brings us to the conclusion that the modeling is not 
sufficient to predict a positron fraction which can be in agreement with AMS data at low 
and high energies at the same time. 

Starting from this observation, in the next Sections, we will extend our model to in¬ 
vestigate whether, and under which conditions, adding an additional positron component to 
the astro model can improve the interpretation of AMS-02 data. Our investigation will be 
oriented in two directions: a contribution coming from dark matter annihilation; a less exotic 
astrophysical interpretations in terms of additional pulsars. 

Before moving on to the discussion, we remark that all the results that are shown in the 
following Sections are obtained by htting only the positron fraction and the total e"*“ -|- e~ 
flux. The reason for this choice is that, as it has been discussed in the previous paragraphs, 
the issue in reproducing AMS-02 regards solely the positron fraction. As we have commented 
in [28] the four observables are the result of independent experimental analyses and therefore 
could be considered to be, at least to some extent, independent. Nevertheless, since under 
a theoretical point of view only two of them are really independent, by htting all the four 
datasets, a certain degree of redundancy enters in the analysis and could hide the information 
concerning the positron fraction. This can be clearly seen in the results that we report in 
Table 1, where the ht with four datasets is associated to a low chi-square despite the tension 
of the model with the measured positron fraction. 

The impact that the choice of the datasets included in the ht has on the results of some of 
the analyses that will be discussed in the following is discussed in Appendix A. 

4.2 Adding dark matter to the picture 

We ht the positron fraction and the total hux within a model that consists of all the primary 
and secondary astrophysical contributions discussed in the previous paragraph plus huxes 
produced by DM pair annihilation or decay. For dehniteness, we refer to this modeling as the 
astro+DM model. The free parameters are now 8: the same Qo.SNRs, 7snrs, AVeia, ^PWNej 
TPWNej (t> already introduced before, plus the two parameters characterizing the DM particle, 
namely its mass m£)M and annihilation cross-section {av) (or decay rate Tdec = I/'T")- 
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Figure 5. Top left panel: the total (statistical plus systematics) uncertainties divided by the central 
value of the data for the four observables measured by AMS-02, are shown as functions of the energy. 
Top right panel: the AMS-02 positron fraction data are reported together with predictions of the 
astro model for different values of the ypwNe parameter, as reported in the boxed inset. The curves 
are normalized in order to be equal at 41.66 GeV. Bottom panel: the associated with the positron 
fraction (obtained in the fit of the positron fraction and total flux datasets) is shown, as a function 
of the energy, for different values of the ?7pwNe and ypwNe parameters. 


Our investigation is twofold, being aimed at: i) deriving bounds on the DM annihi¬ 
lation/decay rate as a function of the DM mass (we will perform this analysis in a raster 
scan of the DM mass); ii) determining whether AMS-02 data require/prefer/allow a DM 
contribution. In both cases, our underlying assumption is that the DM contribution to the 
electron and positron fluxes occurs in the background of emission from astrophysical sources 
(SNRs and PWNe): this is likely the more plausible situation, and we wish to investigate 
whether, to what extent and how robustly a DM leptonic emission can be present in the CR 
fluxes. Our analysis will not restrict the DM contribution to be neither the dominant one, or 
the subdominant one a priori: we will just leave free the contribution from all sources (DM, 
SNRs and PWNe) and adapt the modeling to the data with the same MCMC procedure 
discussed above. 
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Figure 6. Constraints on the DM annihilation cross section {av) (left panel) and DM lifetime r 
(right panel), for various production channels (as reported in the insets) in the astro+DM model. The 
constraints are computed at the 2cr C.L., and refer to the Med propagation model and Einasto DM 
density prohle in our Galaxy. 

4.2.1 Constraints on dark matter annihilation/decay 

Let us move to the determination of upper bounds on the DM annihilation cross section {av) 
(or lower bounds on the decay lifetime r) within a realistic full model for the astrophysical 
emission of primary and secondary e^. We proceed as follows: for every DM annihila¬ 
tion/decay channel and for each value of the DM mass, we perform a MCMC sampling of 
the space defined by the 7 free parameters that characterize our model (6 astrophysical pa¬ 
rameters plus {av) or r). By marginalizing over the astrophysical-sources parameters, we 
derive the posterior distribution functions for the DM annihilation cross section (or lifetime, 
in case of decaying DM), from which 2fT C.L. upper limits are obtained. This refers to the 
values (iTt')iim and rum for which P{{av) < (<Tu)iim) = 97.73% and P{t > rnm) = 97.73%, 
being P the posterior distribution function of the DM annihilation cross section or lifetime. 
The results for both annihilating and decaying dark matter are shown in Fig. 6. 

The bounds are quite constraining, especially in the low DM-mass region. For the 
annihilating case, the softer hadronic channels (here represented by the bb channel) exclude 
thermal annihilation cross section for DM masses below 40 GeV. The same occurs for the 
harder and channels, while for e"*'e“ channel, which has the most prominent 

spectral feature [68] , the bound on thermal cross sections increases at 200 GeV. 

For the leptonic annihilation channels, we notice two features. The first is that for 
lighter DM, the bounds can be quite strong. Annihilation into leptons, especially the direct 
annihilation into a e'’'e“ and to some extent also have a hard energy spectrum with a 

characteristic feature [68], with a sharp fall-off; when the DM mass is below about 100 GeV, 
the feature occurs in the lower-energy portion of the positron flux, where the astrophysical 
contribution is dominated by the soft and smooth secondary flux which, as already mentioned, 
reproduces remarkably well the AMS-02 data. In this case, significantly strong bounds are 
obtained. The limits obtained here are on average in good agreement with the ones that have 
been derived in Refs. [24, 25], even though the astrophysical modeling in our analysis and in 
Refs. [24, 25] differ. 

The second feature observable in the bounds of Fig. 6 for the e^e“ and channels 
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Table 2. Best-fit configurations for the astro+DM model. 


is represented by a series of bumps that appear in the constraints. This corresponds to the 
fact that for some DM mass values, a contribution from DM annihilation or decay is actually 
preferred by the fit, and therefore the bounds weaken. These facts are investigated in the 
next Section. 

4.2.2 Allowed regions in the DM parameters space 

We therefore explicitly investigate which are the configurations in the DM parameters space 
that can provide the best fits to AMS-02 data. Instead of a raster scan on the DM mass, 
we investigate the full 8-fold parameter space of the astro + DM model in order to derive 
the best fit configurations together with their allowed regions. The results are reported in 
Table 2, while the triangular plot of Fig. 7 illustrates the posterior probability distributions 
of the parameters for the case in which DM annihilates into The la and 2cr contour 

plots in the {muM, plane for the annihilating DM case are reported in the upper row 
of Fig. 8 together with the bounds that have been found in Refs. [69, 70] as a result of an 
investigation of the different components contributing to the isotropic gamma-ray background 
(IGRB) measured by the Fermi-LAT experiment. Lastly, in the lower row of Fig. 8 the 
different contribution to the positron fraction measured by AMS-02 are shown for the best- 
fit configurations of the two cases of DM particles annihilating into the and bb channels. 

As it can be clearly seen, the fitting procedure can reach a very high resolution in 
determining the mass and the annihilation cross section/lifetime for the leptonic annihila¬ 
tion/decay channels. This is a consequence of the fact that these channels are expected to 
produce sharply peaked positron spectra which basically contribute to replenish the gap at 
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Figure 7. Triangular plot for the fit of the parameters of the astro+DM model to AMS-02 data, for 
a DM annihilating into the channel. The contours refer to Icr, 2cr and Scr C.L. allowed regions. 

The plots along the diagonal show the posterior distribution for each parameter. 


intermediate energies {i.e. between 20 and 40 GeV), where the secondary flux begins to 
become negligible and the contribution from PWNe starts getting relevant. In this relatively 
small energy window, small variations in DM parameters can lead to large fluctuations of 
the chi-square and this makes the contour regions in the {rriDM, (o’u)) plane quite small. 
A different situation occnrs for hadronic channels, for which the DM contribution extends 
over several energy bins, from intermediate energies up to the last AMS-02 bins: in a large 
fraction of this energy range, the contribution from DM is degenerate with the contribution 
from PWNe and thus quite large variations in DM properties do not always translate into 
large fluctuations on the chi-square. 
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e*/(e'"+e') for DMDM -> nV 


eV(e*+e') for DMDM -> bb 




Figure 8. Upper row: allowed regions in the {judm, {^v)) for a DM particle annihilating into leptonic 
{left panel) or hadronic (right panel) channels. For every channel, the different shadings correspond to 
the la and 2a contour regions, while the solid and dashed lines represent, respectively, the conservative 
and optimistic gamma-ray upper limits derived in Refs. [69, 70]. In the lower row, the contributions 
to the positron fraction associated to the best-fit configurations for a DM annihilating in pf'pT (left 
panel) and bb (right panel) are shown. 


Concerning the comparison with gamma-ray constraints on the same annihilation chan¬ 
nels, the upper row of Fig. 8 reports the upper limits on the DM thermally averaged annihila¬ 
tion cross section (av) computed in Refs. [69, 70] by means of a detail modelling of the various 
contribution to the IGRB energy spectrum that can be derived from the measurements per¬ 
formed by the Fermi-LAT experiment. In particular, solid lines represent conservative upper 
limits as reported in Ref. [69] with the y-ray emission from extragalactic sources set to the 
minimal level allowed in their modeling, while dashed lines refer to optimistic upper limits 
derived in Ref. [70] from a statistical ht to the Fermi-LAT IGRB data. As it can be seen, 
only the annihilation in the e+e“ and pAp,~ channels is fully consistent with the gamma-ray 
bounds, while for all the other channels a tension is present, even in a conservative setup. 

It is worth to emphasize that the astro + DM model has a better agreement with AMS-02 
data than the pure astro model. This is especially true for p~ and r+r” annihilation/decay 
channels, for which the reduced drops to values around 1 (for the astro model we had 
X^/d.o.f = 1.35). In both cases, the predicted positron fraction is in good agreement with 
AMS-02 data: as an example, for the annihilating DM case, the x^ per data point of this 


- 15 - 


















Figure 9. Left panel: positron fraction for the best-fit configurations of the astro model and of the 
astro + DM model p.~ channel), together with AMS-02 data, as a function of the energy. Right 
panel: the associated to the two models is shown as a function of the energy. 


observable is 1.00 for the channel and 1.01 for the channel. To give a better 

illustration of this agreement, we plot in Fig. 9 both the positron fraction for the astro + 
DM model compared with the one predicted by the astro model and the for each data 
point for the two models. As it can be clearly seen, the addition of the DM contribution 
brings a significant improvement of the ht that is not confined to the low-energy data points, 
but extends over all the energy range of AMS-02 measurements. 

Finally, we notice that also for the astro+DM model, the positron fraction keeps a 
flat behaviour beyond the current AMS-02 data: this is a consequence of the fact that the 
addition of a DM component is driven to contribute to intermediate energies, and leaves the 
high-energy part of the positron spectrum basically unchanged from the one obtained in the 
pure astrophysical interpretation. 

4.2.3 Varying the astrophysical modeling 

In this Section we explore the consequences of relaxing some of the assumptions that, in the 
previous Sections, have been made about the contribution from PWNe. This should allow us 
to understand to what extent the remarkable agreement that we have found in the astro+DM 
model can be interpreted as a genuine hint of DM, capable of surviving even within a more 
complex modeling of the PWNe contribution. 

As a first analysis, we carry out a simple exercise to assess what would be the properties 
of a putative pulsar that, added to the ones of the ATNF catalogue, would be able to provide 
an agreement to AMS-02 data as good as the one obtained with the addition of DM. We are 
therefore trying to replace the DM electron-positron production with a generic additional 
PWN emission, to see if this can fully mimic a DM contributions and what are the properties 
required for such a pulsar. To this aim, we perform a fit to AMS-02 data within a version of 
the astro model, modified in such a way to include an additional PWN, whose distance dpsr, 
age Tpsr and total energy emitted in the form of pairs r7psrWo,psr are to be determined by 
a fit to the data. This brings the total number of free parameters in our fit to 9: (5 o,SNRs) 
TSNRs) A^Yela, ^PWNej TPWNej C^psD Tpsj. and ?ypsrkFo,psr- 
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do 
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00 

0 
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0.63l[]:i 

Tpsr [kyr] 

11101873 
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x‘^/84 d.o.f. = 0.85 


Table 3. Best fit for the astro model with an additional PWN. For every parameter the best-ht value 
is reported together with the upper and lower la uncertainty. The value of rjpsr has been derived 
assuming a spin-down luminosity for the additional PWN equal to E = 10^^ erg/s, for definiteness. 




Figure 10. Left panel: the positron fraction for the best-ht conhgurations of the astro model and of 
the astro + additional PWN model are shown together with AMS-02 data as a function of the energy. 
Right panel: the associated to the two models is shown as a function of the energy. 


The result of the fit is shown in Table 3. The left panel of Fig. 10 shows the predicted 
positron fraction for the best-fit configuration compared to the prediction given by the astro 
model, together with the AMS-02 data, while the right panel illustrates the behaviour of the 
as a function of the energy for the two models. The “fictitious” pulsar that we are adding 
sizably improves the agreement with the AMS-02 data as compared to the astro model and 
provides a ht that is comparable to the ones that are associated to a DM signal in the muon 
channel, as reported in Table 2. For the positron fraction dataset, the agreement is extremely 
good: the chi-square per data point is 0.93. This occurs if this source is able to provide a 
non-negligible contribution in the low/intermediate-energy region of the positron flux, which, 
as already stressed, is where also the DM contribution from leptonic annihilation/decay 
channels was preferred by the ht. The best-ht conhguration and the la and 2cj C.L. regions 
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Figure 11. Contour plot of the distance and age of the additional pulsar together with the other 
sources from the ATNF catalogue; the star denotes the best-fit, the darker and lighter shadings 
represent, respectively, the la and 2a confidence regions, while the bars quantify the Icr uncertainties 
on the position of some of the ATNF catalog pulsars (these sources are here circled in red, for clarity). 


for the additional pulsar in the (dpsr, Tpsr) plane are reported, together with the PWNe listed 
in the ATNF catalogue, in Fig. 11. The red bars denote the 1 a uncertainty associated 
to the position of some of the sources that are among the closest to the contour regions 
of the model (as reported in [71]). As one can see, in order to have a seizable flux at 
these energies, the pulsar should be relatively close (600 pc), young (1000 kyr) and powerful 
(??psrkFo = l-4fto4 X 10^® erg). Moreover, the solution appears to fall in a region where the 
ATNF catalogue has no sources. However, once that the uncertainty on the determination 
of the distance of the pulsars is taken into account, some of the sources of the catalogue 
may overlap with the contour regions of the best-fit that we find. In addition, one must 
always remember that this model could be refined by allowing for a free spectral index of the 
fictitious source. This could impact the reconstructed values of the age and distance of the 
putative source. With all this considered, we cannot exclude that this model could represent 
a perfectly viable interpretation of AMS-02 data. 


As a second extension of the analysis of the previous Sections, and in order to investigate 
at a deeper level what we can do with the sources that are in the ATNF catalogue, we 
build a more complex astrophysical framework, which we label refined astro model, that is 
characterized by the assumption that the whole set of PWNe listed in the ATNF catalogue do 
not share anymore the same efficiency r/pwNe- Instead, we single out the five most powerful 
sources of the catalogue in the different energy ranges by using the same ranking algorithm 
that has been used in Ref. [28] and we associate to their emissions the efficiencies r/j (being 
i = 1,..., 5) that act as free parameters in our fitting procedure. These most powerful sources 
are [28]: Geminga (i = 1), J2043-I-2740 (i = 2), J0538-I-2817 {i = 3), Monogem {i = 4) and 
B1742-30 (z = 5). To limit the number of parameters in the fit we associate to each one of 
these five sources the same spectral index 7 i=i ,„.,5 = 1.8 [28] (we have checked that different 
values for this parameter does not affect our conclusions). The parameters ?7pwNe and ypwNej 
whose definition has been given in the previous Sections, are still present and characterize 
the population of the other PWNe in the ATNF catalogues, which behave like a background 
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Table 4. Best-fit configurations for the refined astro model. 


of sources that share common values for their efficiencies and spectral indices. To summarize: 
we consider emission from the full ATNF catalogue, but for the 5 most powerful sources we 
allow for a free normalization parameter on their emission. 

In Table 4 we list the best-fit configurations that are obtained by fitting AMS-02 data 
within the refined astro model, without and with the addition of DM. We illustrate only the 
results obtained for the annihilating case, since we do not register any significant change for 
a DM that decays in the same channels (except for the fact that the best-fit DM masses are 
halved). We show only the best-fit configurations for the leptonic channels, since for DM 
annihilation into the bb and W~^W~ channels the addition of a DM signal does not improve 
the fit as compared to the purely astrophysical case. 

In the case of a pure astrophysical interpretation, the results in terms of the recon¬ 
structed parameters are stable as compared to the ones obtained in the astro model, with the 
notable exception of the emission coefficient r]i for the 5 most powerful PWN: Geminga and 
Monogem are required to be the dominant emitters. The catalogue-PWN other than the 5 
most powerful are now significantly suppressed as compared to the previous analysis of the 
astro model. riPWNe drops from 3.7% to 0.1%. This implies that the AMS-02 data could be 
shaped by very few dominant sources, like Monogem and Geminga. The significance of the 
fit is here improved as compared to the simpler astro model and makes this model perfectly 
compatible with AMS-02 data. This holds also for the positron fraction dataset for which 
the chi-square per data point is 1.02. 

In the case a DM signal is present, we see that within this more complex model for 
the astrophysical emission, the contribution from a DM annihilating into leptons is able to 
improve only very slightly the fit to AMS-02. The chi-square per data point of the positron 
fraction is reduced only for the annihilation channel, and reaches 0.9. Concerning 

gamma-ray upper limits, the best-fit configuration for the case is still in tension with 

the constraints shown in Fig. 8, while the DM candidate is fully compatible with these 

bounds, and therefore represent the best option for a DM interpretation. In the case of the 
preferred /ichannel, the best-fit mass in the refined astro model turns out to be almost 
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Figure 12. la and 2a allowed regions in the {itidm, (o'v)) plane for a DM annihilating in the 
yT channel in the refined astro model with DM. Solid and dashed lines represent, respectively, the 
conservative and optimistic gamma-ray upper limits derived in Refs. [69, 70]. 


eh(e*+e') POWERFUL PWN + DMDM -> nV 



E [GeV] 


Figure 13. The contribution to the positron fraction associated to the best-fit configuration for a 
DM annihilating in channel (solid green line) shown together with the flux from the five most 

powerful sources (dot-dashed magenta line), from all the other pulsars of the ATNF catalog (dotted 
blue line) and from secondary production (dashed red line) 


the same as in the astro model (82 GeV instead of 89 GeV) and the annihilation cross section 
is almost a factor of two smaller, being compatible with the thermal value at the 2cr level. 
The allowed region and the gamma-rays bounds for the channel are shown in Fig. 12. 

We notice also that in this case, Monogem is the largely dominant pulsar emitter. 

The positron fraction corresponding to the best fit results for the refined astro model 
with DM, which occurs for the channel, is shown in Fig. 13. Notice that in this case, 

where a single pulsar largely dominates the positron flux (Monogem), the positron fraction 
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is predicted to drop relatively quickly at the pulsars cut-off energy (which in our analysis is 
taken at 2 TeV). This feature is here due not to the DM emission (which, in this case, gives 
its maximal contribution at intermediate energies), but to the pulsars high energy properties. 
Data above the TeV scale could help in distinguishing the situation where all known pulsars 
contribute in a democratic way {astro model, with or without DM), from a case where a 
single pulsar dominates the high-energy positron emission (in our analysis, the refined astro 
model, again with or without a DM contribution). 

5 Conclusions 

In this paper we have performed a quantitative study of AMS-02 electron and positron data, 
in terms of both galactic known sources and a possible dark matter contribution. 

Firstly, by updating the analysis illustrated in Ref. [28], we have investigated the pos¬ 
sibility to interpret the whole set of leptonic AMS-02 data in terms of astrophysical sources 
of primary and secondary origin only. We have shown that a model based on the simplifying 
assumptions that all the PWNe of the ATNF catalogue share the same efficiency and spectral 
index is in some tension with the measured positron fraction. 

Secondly, we have carried out a detailed study of the interplay between the contributions 
to the positron flux that derive from primary astrophysical sources (namely PWNe) and 
from the annihilation/decay of DM particles. We have worked within a scenario in which 
DM gives its contribution on top of a realistic astrophysical leptonic emission, which acts 
as a non trivial background with respect to a putative DM signal. On one hand, we have 
investigated how one can use the highly precise data provided by the AMS-02 experiment 
to derive robust constraints to DM properties. On the other hand, we have illustrated 
how the addition of a DM contribution compatible with the bounds that arise from other 
indirect detection channels, can improve the fit to AMS-02 data as compared to the simple 
astro model. Specifically, we have found that a DM particle with a mass around 80 GeV, 
annihilating in the channel with a cross section remarkably close to the thermal value 

can provide, when added to the PWNe and to the secondary contributions, an excellent fit to 
the AMS-02 data, positron fraction included. Conversely, the addition of a DM annihilating 
or decaying into hadronic channels, does not appear to improve the fit significantly and 
the best-fit configurations for these channels are in tension with upper bounds derived from 
different indirect DM searches (e.g. gamma-rays). 

We have also investigated whether the addition of a local pulsar to the sources of the 
ATNF catalogue can provide a fit to the AMS-02 data of comparable goodness as the DM 
hypothesis. This is indeed possible and, within this theoretical framework, one is able to 
obtain a remarkably good fit to AMS data. The characteristics (distance and age) of this 
additional pulsar do not correspond exactly to any source already present in the ATNF 
catalogue. However, once that the uncertainty associated to the distance of these sources is 
taken into account, some of them appear to possess characteristics that are compatible with 
the best-fit region of the model. This makes this interpretation of AMS-02 data perfectly 
viable. 

Finally, we have shown that relaxing the hypothesis of equal efficiency and spectral 
index for all the PWNe of the ATNF catalogue weakens the need to have additional sources, 
and a pure astrophysical interpretation of AMS-02 data is indeed possible. 

We conclnde by saying that the parallel inspection of independent and multi-wavelength 
channels is a crucial step in pursuing a deeper understanding of the astrophysical sources that 
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populate our Galactic environment. This investigation, together with the precise measure¬ 
ments performed by current and future experiments may significantly help in interpreting 
the cosmic lepton data and, perhaps, in shedding light on the DM mystery. 
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A Fitting different datasets 

As discussed in Section 4, all the results shown in this paper have been derived by fitting 
only the positron fraction and total flux datasets. In fact, as we have already discussed, since 
the four observables measured by AMS-02 are not fully independent, including additional 
datasets in the fit could potentially lead to an overfit of the data and bias our results. 

To illustrate this point, we plot in Fig. 14 the la and 2a contour regions for the best-fit 
parameters of different models obtained by fitting different sets of data (we show the astro + 
DM model in the top row and the astro model with an additional PWN in the bottom row). 
In particular, the blue contours are the ones obtained by fitting the four datasets, while the 
red ones are the result of fitting only the positron fraction and the total flux. It is manifest 
that fitting all the four datasets determines a shrinking of the contour regions. The best-fit 
values of the parameters of the different models that we consider seem to be not significantly 
affected by the choice of the datasets to include in the fit, but the value itself of the chi- 
square changes notably. In fact, when fitting the four datasets, the best fit configurations are 
associated to yfl/d.o.i. = 0.78 (astro + DM model, pflpT channel), y^/d.o.f. = 0.84 (astro + 
DM model, r+T“ channel), y^/d.o.f. = 0.70 (astro model -|- additional PWN) while, for the 
same models, we have y^/d.o.f. = 0.98, y^/d.o.f. = 1.05 and y^/d.o.f. = 0.85 when fitting 
only the positron fraction and the total flux. Usually, values of the reduced chi-square that 
are significantly less than 1 represent a clear indication that the choice of datasets is such 
that the model is overfitting data. This is the reason why, for all the analyses shown in the 
paper, we have considered fits to the positron fraction and the total flux only. 

Lastly, in Fig. 15, we show how the choice of the datasets to be included in the fit affects 
the determination of the upper limit to the DM annihilation cross section for the 
annihilation channel. As one could have easily expected, in connection with the shrinking 
of the contour regions discussed in the previous paragraph, for the parameters around the 
best-fit ones, the bounds obtained with fits to only the positron fraction and the total flux 
are weaker than the ones obtained by fitting the four datasets. On the other hand, for all 
the masses for which the addition of DM does not improve the quality of the fit with respect 
to the astro model, bounds are rather insensitive to the choice of dataset. 
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Figure 14. la and 2a contour regions for different models and different choices of the datasets 
included in the fit. Blue lines refer to the fit of all the four AMS-02 datasets, while red lines refer 
to the fit to the positron fraction and total flux only. Stars denote the best-fit configurations. The 
contour regions reported in the top row refer to the astro -I- DM model for two different annihilation 
channels: in the left panel and r+r” in the right panel. The bottom row shows the contours 

in the (d,T) plane for the additional PWN discussed in Section 4.2.3. 



Figure 15. Upper limits to the DM annihilation cross section, for the annihilation channel, 

for different choices of datasets included in the fit, as reported in the boxed inset. 
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